#include <iostream>
#include <unsupported/Eigen/MatrixFunctions>

using namespace Eigen;

int
main()
{
	using std::sqrt;
	MatrixXd A(3, 3);
	A << 0.5 * sqrt(2), -0.5 * sqrt(2), 0, 0.5 * sqrt(2), 0.5 * sqrt(2), 0, 0, 0, 1;
	std::cout << "The matrix A is:\n" << A << "\n\n";
	std::cout << "The matrix logarithm of A is:\n" << A.log() << "\n";
}
